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ABSTRACT 

Context. A discharge of nitrogen gas, as created in a microwave-induced plasma, exhibits a very dense molecular emission line 
spectrum. Emission spectra of this kind could serve as wavelength calibrators for high-resolution astrophysical spectrographs in the 
near-infrared, where only very few calibration sources are currently available. 

Aims. The compilation of a spectral line list and the characterization of line intensities and line density belong to the initial steps when 
investigating the feasibility of potential wavelength calibration sources. Although the molecular nitrogen spectrum was extensively 
studied in the past, to our knowledge, no line list exists that covers a continuous range of several thousand wavenumbers in the 
near-infrared. 

Methods. We recorded three high-resolution (Ay = 0.018 cm"^) spectra of a nitrogen gas discharge operated at different microwave 
powers. The nitrogen gas is kept inside a sealed glass cell at a pressure of 2 mbar. The emission lines in the spectra were fitted by a 
superposition of Gaussian profiles to determine their position, relative intensity, and width. The line parameters were corrected for an 
absolute wavelength scale, instrumental line broadening, and intensity modulation. Molecular and atomic transitions of nitrogen were 
identified with available line positions from the literature. 

Results. We report line lists with more than 40 000 emission lines in the spectral range 4500-11000 cm"^ (0.9-2.2 pm) *. The spectra 
exhibit emission lines over the complete spectral range under investigation with about 350 - 1300 lines per 100cm"k Depending on 
the microwave power, a fraction of 35% - 55% of all lines are blended. The total dynamic range of the detected lines covers about 
four orders of magnitude. 

Conclusions. Line density and peak intensities qualify the recorded discharge as a useful wavelength calibrator, and the line list 
provides an empirical reference for nitrogen spectra in the near-infrared. 

Key words. Techniques: spectroscopic - Methods: laboratory: molecular - Methods: data analysis 


1. Introduction 


High-resolution near-infrared spectroscopy is an increasingly 
important tool in astrophysics. A major driver for this devel¬ 
opment is research on cool low-mass stars, such as M dwarfs, 
which emit most of their photons at wavelengths above 1 pm. 
Eor example, M dwarfs are of great importance in the search for 
low-mass extrasol ar planets and are in the focus of upcoming ob - 
serving programs ( [Reiners et al.|2010]|Quirrenbach et al.|2014| ). 
A key element of precision spectroscopy in astronomy is a reli¬ 
able wavelength calibration, which assigns the wavelength scale 
to the pixels of the spectrograph’s detector. 

Currently, the two most frequently used wavelength calibra¬ 
tion techniques are absorption gas cells (e.g., filled with iodine) 
or hollow cathode lamps (e.g., [Marcy & Butl^|1992t [Baranne 
et al.|1996| ). Emission lamps with atomic line spectra have been 
used at optical wavelengths with great sucess, but they feature 
fewer lines toward longer wavelengths. The line list for commer¬ 
cial thorium-argon hollow cathode lamps contains about 2400 
spectral lines suitable for wavelength calibration in the range 
900 - 4500nm ( jKerber et al.||2008| ). Absorption cells with im¬ 
proved gas mixtures, and new techniques, such as laser fre¬ 


* The line lists and spectra are available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsarc.u-strasbg.fr/viz-bin/qcat7J/A-hA/vol/page 


quency combs and Eabry-Perot etalons, are under development 
(jSeemann et al.jjMT^ jWilken et al.||2012l [Schafer & Reinersj 
[2012[ ). Nevertheless, molecular emission spectra, e.g., from N 2 
or CN, can be a viable alternative for wavelength calibration 
( Boesch et al.|2014| ): a discharge of nitrogen gas exhibits a very 
dense line spectrum in the near-infrared, in contrast to hollow 
cathode lamps; the frequency of all lines are defined by physical 
laws not supposed to change over time, in contrast to Eabry-Perot 
etalons; and the equipment is relatively cheap, in contrast to laser 
frequency combs. 

Eor the investigation of new calibration sources, the compi¬ 
lation of a spectral line list and the characterization of line in¬ 
tensities and distribution are initial steps. Although the spectrum 
of molecular nitrogen has been extensively studied in the past, a 
line list based on high-resolution observations over a wide range 
in the near-infrared, to our knowledge, does not exist. 


An extensive compilation of spectroscopic data on the 
molecule N 2 (and its ions N“, NJ and N^^) was presented by 
Lofthus & Krupeni^ ( [1977[ ). Their review discusses each elec¬ 


tronic band system, provides extensive references to previous re¬ 
search, and can be used to identify band systems expected in a 
given wavelength range. Eor N 2 , the Eirst Positive System B^Hg- 
A^E^ (478 - 2531 nm) is the most prominent band system and 
appea r s in m ost types of discharges according to [Lofthus & Kru- 


penie 


(1977). Eurthermore, the Meinel System A^Hu-A^Sg of 
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NJ has transitions in the near-infrared (550 - 1770 nm). More 
recent reports on spectroscopic data for molecular nitrogen in 
general can, e.g., be found in Laher & Gilmore| (l 991 1 ) , [GTlmore 
|et al. ( [1992| ) and Laux & Kruger] ( | 199^ 


The First Positive System was specifically analyze d by |Ef- 
fantin et ar](|1979| ) and |Roux et al.| ( |198^ . The paper by |Effaiitin 


et al.| ( |197'^ covers the vibrational bands (0-0), (1-0) and (2-0) 


and includes line lists of observed and calculated wavenumbers. 
The study of the First Positive System was later extended to 33 
bands by |Roux et al.| ( [1983] ), but they did not publish their line 
list. Our efforts to retrieve the list from the journal or the insti¬ 
tute where the research was conducted were not successful. For 
the Meinel System of , line lists of observed and calculated 
wavenumbers for the (0-0), (0-1), (1-0) and (1-2) bands are in¬ 
cluded in the work by |Ferguson et al.| ( |l992 ). 

In this article, we present the characterization of a nitrogen 
gas discharge observed in the near-infrared with a Fourier trans¬ 
form spectrometer (FTS). Three spectra were recorded with the 
discharge operated at different microwave (MW) powers. We de¬ 
scribe the experimental procedure and the data reduction pro¬ 
cess in detail. We fit all emission lines above a minimum flux 
threshold and produce a list including line position, line inten¬ 
sity, and line width. The identification of all transitions based on 
theoretical calculations is beyond the scope of this work, but we 
include assignments of molecular nitrogen transitions from the 
(0-0) band of the First Positive System and from four bands of 
the Meinel System, in addition to atomic nitrogen lines. We an¬ 
alyze properties of line intensity and line density with a focus 
on the spectrum’s usefulness as a wavelength calibrator in as- 
trophysical spectroscopy. For the discharge operated at different 
MW powers, we investigate the change in line width of detected 
lines and in line intensity of identified atomic nitrogen. 


2. Experiment 


The experim ental setup is similar to that described in |Boesch| 
et al.lpOT^ . Fig. [2 shows a schematic drawing. Here, we sum¬ 
marize the general functionality of the experiment for complete¬ 
ness and highlight the details that are important for the measure¬ 
ments presented here. 

We use an Evenson-type microwave cavity (from Opthos In- 
strument0 and a solid-state microwave generator (GMS 200W 
from Sairem) to create a glowing plasma inside a sealed dis¬ 
charge gas cell. The cell has a diameter of 12 mm, a length of 
about 20 cm, and is made of quartz glass. A plane quartz win¬ 
dow is melted to the front end with a tilt of 5°. The cell, filled 
with 2 mbar of high-purity dry nitrogen gas, was purchased from 
Sacher Lasertechnik. A sealed gas cell enables us to make re¬ 
peated measurements with the same gas properties and allows 
for a compact experimental design. 

The emitted light passes through a convex lens (Thorlabs 
LA1255) and is focused by a reflective collimator (Thorlabs 
RC08SMA-P01) onto an optical fiber (Thorlabs EG550LEC). 
The fiber feeds the lights into an ETS (Bruker lES 125HR). On 
Eebruary 14, 2015, we recorded three high-resolution spectra 
(Ay = 0.018 cm“^; spectral resolving power ofR = 555556 at 
A = 1 pm) of the nitrogen gas discharge. Eor each spectrum, 
the discharge was operated at a different MW power. The for¬ 
ward MW power set at the generator was 25 W, 50 W, or 100 W 

^ Commercial products are identified in this document to specify the 
experimental procedure adequately. This identification is not intended 
to imply recommendation or endorsement, nor is it intended to imply 
that the products identified are necessarily the best available for the pur¬ 
pose. 
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Fig. 1. The experimental setup. See main text for explanation (Sect.[^. 


Table 1. Settings of the FTS for each of the three measurements with 
different MW powers. 


resolution 

0.018 cm-' (x0.944‘') 

scans 

300 (j^ 6 h 45 min total scanning time) 

aperture 

1.3 mm 

detector 

InSb (liquid nitrogen cooled) 

beamsplitter 

CaF2 

pressure inside FTS 

0.10 - 0.15 mbar 

zerofilling factor 

4 


Notes. The applied apodization (Norton-Beer medium) modifies the 
resolution set in the spectrometer control software by this factor. 


for the three measurements. The chosen resolution allowed us to 
record a symmetric interferogram that minimizes phase correc¬ 
tions errors. For each of the three measurement runs, the final 
interferogram is an average of 300 scans. The parameters of the 
FTS are given in Table[2 

For the determination of the absolute wavenumber scale, we 
additionally recorded spectra of a sealed discharge cell filled 
with 2 mbar of argon gas. The argon gas cell simply replaced 
the cell filled with nitrogen for these measurements. 

Flat-field spectra were measured using a laboratory halogen 
lamp (Ocean Optics HL-2000-FHSA-HP) to correct spectral line 
intensities for the transmittance of optical elements and the de¬ 
tector response. For the fiat-field measurements, the gas cell was 
removed from the experimental setup and the halogen lamp was 
placed behind the MW cavity. 

Figure shows the raw spectrum of the nitrogen gas cell, 
covering the wavenumber range 4000 - llOOOcm"^ (0.9 - 
2.5 pm). Intensities are given in arbitrary units (a.u.) as obtained 
from the Fourier transform. The transmittance of the optical fiber 
and the response curve of the detector set the limits of the usable 
spectral range. No signal above noise was detected at wavenum¬ 
bers below 4200 cm“^ (above 2.38 pm). This part of the spectrum 
is indicated in Fig. with a shadowed box and was used for de¬ 
termining the absolute noise level (see Sect. . The spectra are 
available in the online material. 


3. Data Analysis 

This section explains the steps in the data reduction to obtain the 
parameters of the spectral emission lines. First, we determine the 
absolute noise level in the spectrum. Then, the emission lines are 
located and fitted. The line parameters are corrected for absolute 
wavenumber scale, instrumental broadening, and intensity mod- 
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wavelength [pm] 

2.5 2 1 0.9 



Fig. 2. Raw emission spectrum of the nitrogen 
gas discharge operated at a MW power of 50 W. 
Strong vibrational bands from the First Positive 
System of are labeled (v' - v") for orien¬ 
tation using band origin wavenumbers as given 
by |Gilmore et al.|(T992| ). The spectral range be- 
low 4200 cm" ^ (gray box) does not contain any 
emission lines and was used to determine the 
absolute noise level in the spectrum. 


ulation. Finally, spectral lines are identified with available line 
positions from the literature. 


3.1. Noise level 


The knowledge of the noise in the spectrum is important to prop¬ 
erly weight the data points and derive meaningful errors for the 
line parameters through error propagation. Under the assumption 
of white noise in the interferogram, the noise in a spectrum ob 
tained from an FTS is uniformly distributed ( [Davis et al.|2001 


IS 


jVoigtman & Winefordner||1987| ). The premise of white noise 
reasonably valid for the dominance of detector or photon noise. 
The determination of the noise in a single spectrum without 
a sufficient continuum poses a practical difficulty. In our case, 
the spectrum has a baseline without any spectral lines only be¬ 
low 4200 cm“^ (marked with a gray box in Fig.[^. We use this 
wavenumber region to calculate the noise. 


Additionally, the baseline of the spectrum has a small offset 
from zero, which is best visible in the region without spectral 
features. The offset can be explained with an uncertainty in the 
first component of the Fourier transform, i.e., the intensity of 
the interferogram at zero path difference. Because this value is 
wavelength independent, we subtract the median flux from the 
wavenumber region 4000 - 4200 cm“^ from the whole spectrum 
before further analysis. 

Then, we arrange the flux values from the wavenumber re¬ 
gion 4000-4200 cm“^ in a histogram and fit a Gaussian function 
to it describing the distribution very well; Fig. shows this ex¬ 
emplary for the measurement of the discharge operated at a MW 
power of 50 W. We adopt the standard deviation of the Gaussian 
profile as the absolute noise level. The first two rows in Table 
list the subtracted offset values and the noise levels for the three 
measurements. 



Fig. 3. Histogram of the flux values in the spectral region 4000 - 
4200 cm" ^ from the measurement of the nitrogen discharge operated 
at a MW power of 50 W. A Gaussian profile is fitted to the histogram. 

Table 2. Baseline offset, absolute noise level, minimum peak intensity 
threshold, and number of located emission peaks for the three measure¬ 
ments. 



25 W 

50 W 

100 w 

baseline offset [a.u.] 

0.00009 

0.0001 

0.0001 

noise level [a.u.] 

0.0001 

0.0001 

0.0001 

minimum peak intensity [a.u.] 

0.001 

0.001 

0.001 

number of peaks found 

40408 

51776 

58274 


3.2. Peak finding and peak fitting 

We locate the emission peaks in the spectrum, together with the 
intensity at these points, using the zero-crossings in the spec¬ 
trum’s first derivative. The derivative is calculated at every data 

Article number, page 3 ofpfi] 



































































A&A proofs: manuscript no. 26353_final 


point after the entire spectrum has been interpolated with a cubic 
spline. To be classified as a peak, a data point must fulfill the fol¬ 
lowing three criteria: (i) a zero-crossing in the first derivative, (ii) 
an intensity in the spectrum above a threshold value of ten times 
the absolute noise level, and (iii) a maximum intensity within 
one resolution element. Experience shows that the adopted min¬ 
imum intensity threshold value is a good compromise between 
considering as many peaks as possible and satisfactory fit re¬ 
sults. The minimum peak intensity (ten times the absolute noise 
level) and the number of located peaks are given in Table The 
peak search and the final line list cover the wavenumber range 
4500- 11000 cm-i. 

The first derivative’s zero-crossing technique delivers esti¬ 
mates of the lines’ positions and intensities. The estimates serve 
as initial values for the subsequent fitting process, which is 
needed to determine accurate line parameters. The fitting pro¬ 
cedure employs a weighted least-squares minimization (package 
mpfit for the Python programming languag^. Multiple emission 
lines are simultaneously fitted with a superposition of Gaussian 
profiles. The combined fit of several spectral lines is important 
to account for overlapping lines. 

The fit with a Gaussian profile gives the line position vq, line 
intensity /, and line width ctq (standard deviation). We extract 
the uncertainties for the parameters from the covariance matrix 
of the fitting routine. The full width at half maximum (FWHM) 
of a line can be calculated using FWHM = 2o-q ^J2\n(2). We 
also tested a Voigt profile for the fit function, but did not obtain 
better values for the reduced compared to the Gauss profile. 
We rather observed that the Voigt profile is strongly influenced 
by (blended) low-intensity lines near the noise level because the 
overall shape of the Forentzian profile is more sensitive to data 
points in the wings of the lines. 


3.3. Correction of absolute wavenumber scale 

One advantage of Fourier transform spectroscopy is the high rel¬ 
ative precision of the wavenumber scale. The absolute wavenum¬ 
ber scale as given by the internal reference laser is, however, 
subject to changes due to misalignments and environmental 
changes inside the FTS (see, e.g., [Griffiths & de Haseth|2007 


chap. 2.6). The resultant shift in the wavenumber scale is lin¬ 
ear in wavenumber and can be corrected using known positions 
of calibration lines. In theory, one reference line with a well- 
known line position is sufficient, but in practice a set of lines is 
used. The correction factor Kq^ is then the weighted average of 
K = (vo, ref /vo, measured) “ 1 for each reference line, where vq, ref 
is the line position from the literature and vq, measured is the line 
position in the recorded spectrum. The following calibration is 
applied to all detected spectral lines: 


V0,C = Vo(l +/^eff), 


( 1 ) 


where yo,c is the calibrated wavenumber and vq is the uncali¬ 
brated wavenumber. 

To determine the absolute wavenumber scale, we compared 
positions of Ar I li nes with data from the NIST database 
( Kramida et al.|201^ . We recorded spectra of an argon gas dis¬ 
charge before and after the measurements of the nitrogen spec¬ 
tra. The line positions utilized from the NIST database are given 
with an accuracy of three decimal places. Additionally, taking 
information from the original reference of the argon transitions 
( Norlen|[T973| ) into account, we conservatively assume an un¬ 
certainty of 0.003 cm“^ The errorbars on k for the individual 


http://cars.uchicago.edu/software/python/mpfit.html 
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Fig. 4. The correction factor k for the absolute wavenumber scale of 
the FTS as determined for 22 atomic argon lines. The horizontal line 
indicates the weighted average (solid line) with its error margin (dashed 
lines). Shown are the data for the argon spectrum recorded before the 
measurements of the nitrogen spectra. 


argon lines (as shown in Fig. are calculated by error propa¬ 
gation of the uncertainties of the position of the reference line 
and the line’s position in our spectra as determined from a line 
fit. Using 22 argon lines, we calculated correction factors of 
^efF, before = (“0.9 + 0.8) • 10“^ and /Ceff, after = (“0.8 + 0.8) • 10“^ 
from the observations taken before and after the measurements 
of the nitrogen spectra, respectively. Combining the two results, 
we accept a value of /Ceff = (-0.85 + 0.85) • 10“^. This re¬ 
lates to an uncertainty in the corrected wavenumber scale of 
(0.85 • 10“^)yo, namely about 0.00038 cm“^ at 4500 cm“^ and 
about 0.00094 cm“^ at llOOOcm"^ 

3.4. Correction for instrumental resolution 

The shape of a recorded spectral line is the convolution of its 
true shape and the instrumental line spread (IFS) function. We 
assume the IFS function to be of Gaussian shape. The convolu¬ 
tion of two Gaussian functions results in another Gaussian with 
a variance that is equal to the sum of the variances of the origi¬ 
nal functions. Under the assumption of a Gaussian IFS, we can 
estimate the FWHM before the instrumental broadening by ap¬ 
plying 

FWHMc = .y/pWHM^ - FWHMf„^„., (2) 

where FWHM and FWHMc are the uncorrected and corrected 
values, respectively. The FWHMinstr = (0.018-0.944) cm“^ is the 
theoretical FWHM of the IFS function (compare with the reso¬ 
lution given in Table[^. Using FTS observations of a frequency- 
locked diode laser with a line width of 100 kHz, we found 
that the true IFS tends to be slightly broader than the theoreti¬ 
cal value, but the difference is less than 10%. 

3.5. Correction of line intensifies 

The intensity of the spectrum is modulated because of the trans¬ 
mittance of the optical elements in the experimental setup and 
the response curve of the detector. The true intensities of the 
emission lines can be retrieved through comparison with a spec¬ 
trum of known shape, i.e., a fiat-field spectrum. The spectrum 
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Fig. 5. The spectrum of a halogen lamp is used as a flat-field (spectrum 
shown in light blue with absorption bands from atmospheric water). 
The dark blue line describes the continuum of this spectrum. The green 
dashed curve shows the theoretical normalized curve of a blackbody at 
the halogen lamp’s bulb temperature of 3000 K (scaled by a factor of 
1/50 for visualization). The red dot-dashed curve is the flat-field con¬ 
tinuum divided by the theoretical normalized blackbody spectrum. The 
raw intensities of the spectral lines are corrected by the latter curve. The 
vertical dashed line indicates 4500 cm"^ 


of a laboratory halogen lamp was recorded before and after the 
measurements of the nitrogen gas cell, and the two spectra were 
averaged to give one flat-field spectrum (see Fig. 0 . Theoreti¬ 
cally, the halogen lamp produces a blackbody spectrum with a 
maximum intensity at 969 nm, which corresponds to the bulb 
temperature of 3000 K, as referenced by the manufacturer. The 
spectral shape of the blackbody is not constant over the whole 
wavenumber range under investigation (green dashed curve in 
Fig. We correct the fiat-field spectrum’s continuum by the 
theoretical curve and use the result for the intensity correction: 
the line intensities 1 derived from the spectra are divided by the 
intensity of the corrected fiat-field curve at their position (red 
dash-dotted line in Fig.[^. As result, we obtain the fiat-field cor¬ 
rected peak intensities 1^. 

The degree of knowledge of the lamp’s true spectrum and the 
deviation from the assumed theoretical blackbody limit the ac¬ 
curacy in the fiat-fielding process. The spectral output of a sim¬ 
ilar lamp (Ocean Optics HL-2000-FHSA) was measured by the 
manufacturer with respect to a calibrated system at wavenum¬ 
bers above 6250 cm“^ (up to 1.6yum; private communication). 
The maximum deviation of the measured spectral output from a 
theoretical blackbody with the temperature as stated in the spec¬ 
ifications was less than 7%. Taking further uncertainties in the 
fiat-fielding process into account, we assume that the error on 
/c is around 10 % and therefore small compared to the intensity 
modulation by the detector response, which can be inferred from 

Fig-S 

3.6. Atomic line identification 

All emission peaks in the spectrum were checked against line 
lists f rom the atomic line database from NIST ( [Kramida et al.| 
|2014| ) for N I, N II, O I, and Ar I. We checked for oxygen and 
argon because they are the most abundant species in the atmo¬ 
sphere after nitrogen, and would be obvious contaminators, if 


the cell was not filled and sealed properly. We only used lines 
with observed wavelength values in the NIST database. For a 
match, we require the distance between the recorded line and 
the database to be less than the line’s half width at half maxi¬ 
mum (HWHMc = FWHMc/2). The probability that the position 
of a line from the literature matches a line in our recorded spec¬ 
trum by chance is 20%, 30% and 38% for the measurements with 
MW power of 25 W, 50 W and lOOW, respectively, because of 
the many densely-spaced spectral lines. 

Although up to 35 oxygen lines (53% of the lines listed in 
the NIST database) and 35 argon lines (57%) fulfill the above 
mentioned HWHMc criterion, we judge them to be misidentified 
in the crowded spectrum. We come to this assessment, because 
the determined FWHMc of the lines do not follow the expected 
linear trend with wavenumber (see Sect. |4.4| and compare with 
Fig. and the peak intensities Ic do not change by the same 
amount when the MW power is altered (compare with Fig. 12). 
The absence of lines from oxygen and argon indicates that no 
leakage of air into the gas cell had occurred. 

For nitrogen, only lines from the neutral atom could be iden¬ 
tified with confidence. This is not surprising because the degree 
of ionization is low in this kind of electrodeless discharge lamp 
( [Thorne et al.|1999] ). 


3.7. Identification of molecular nitrogen lines 

We use the spectral line lists with branch designations for molec¬ 
ular nitrogen transitions in the wavenumber range of our spectra 
from Effantin et al.| ( |1979] ) and jFerguson et al.| ( 1992] ). As for the 
atomic line identification, we require that the distance to the lit¬ 
erature value is less than the line’s half width at half maximum. 
The work by [Elfantin et al.j ( |1979| ) provides line lists for 


the (0 - 0 ) band of the B'^Hg-A^E^ system of ^'^N 2 and ^^N 2 . 
The lists consist of measured wavenumbers for 21 branches 
and calculated wavenumbers for 27 branches. The observed and 
calculated wavenumbers differ by less than 0.01 cm“\ where 
both values are given. [Effantin et al.j ( [1979[ ) did not observe 
lines from the six branches designated R 13 , Q 13 , P 13 , R 31 , Q 31 , 
and P 31 . Our nitrogen gas contains the natural ratio of isotopes 
(14n2 :1^N2 = 99.636 : 0.364; [Lide[[2Q0^ and therefore only 
very faint lines of ^^N 2 are expected, if any at all. We do not 
include them in the analysis. Eor ^'^N 2 , we can identify up to 
80% of the lines from the literature, including some from the 
above mentioned branches that were not observed by [Effantin 
[eraLl ( [T979l ). 

Eurthermore, we identified up to 89% of the lines from four 
vibrational bands of the Meinel System of ^"^N^ as listed in 
[Eerguson et al.[ ([1992J. The four bands have their origins at 
9016 cm-^ (0- 0), 6841 crn'i (0 - 1), 10889 crn'i (1 - 0), and 
6572 cm-i (1-2). 

The HITRAN2012 database ( [Rothman et al.[2013] ) lists addi¬ 
tional 584 quadrupol transitions in the electronic ground state of 
^^N 2 (in the wavenumber range 4500-11000 cm“^). The number 
of matches in line position (17%, 26%, and 34% for our three 
measurements with MW power of 25 W, 50 W, and lOOW, re¬ 
spectively) is below the probability for a random match for each 
spectrum, as stated in the previous section. 


3.8. Examples of line fitting 

We present in Eig. a part of the three spectra taken at differ¬ 
ent MW powers. Each of the three panels shows the same spec¬ 
tral interval (9892 - 9895 cm“^). The black crosses are the data 
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points of the spectrum and the blue curve represents the fit to 
the emission lines. The black dashed horizontal bar indicates the 
minimum peak intensity (peak threshold). The peaks are marked 
at their line center and maximum line fiux, as determined by the 
fit. A spectral line is marked with a filled symbol if no noticeable 
problems during the analysis process occurred, while an empty 
symbol indicates a line with one or more fiags (see Sect. |4.1| ). 
A green square identifies an N I line and a gray circle marks 
an unidentified line. The lower part of each panel shows the fit 
residuals (blue curve). 

We observe in Fig. that line intensities grow with increas¬ 
ing microwave power. At the same time, the noise level remains 
constant (compare Table [^, and therefore more data points ex¬ 
ceed the peak threshold and are included in the line fit. 


4. Results 

In the following section, we present the spectral line lists created 
from the measurements of the nitrogen gas discharge and some 
properties of the line parameters. The analysis of line intensi¬ 
ties and line density focuses on the usability of this light source 
for wavelength calibration of astrophysical spectrographs. In 
Sect. |4.4| we demonstrate the influence of the MW power level 
on the line width of detected lines in general, and, specifically, 
on the peak intensities of identified N I lines. 

4.1. Spectral line list 

We created separate line lists for the three measurements that in¬ 
clude all detected emission peaks with their parameters in the 
wavenumber range 4500 - llOOOcm"^ (0.9 - 2.2pm). The line 
lists are available in the online material as simple text files con¬ 
sisting of 15 columns. Tablegives a description of the content 
in each column. The stated errors of the positions cr(yo), line 
intensities cr(I), and line widths cr(crQ) are relative errors corre¬ 
sponding to the line fit in the spectrum. They are taken from the 
covariance matrix of the weighted least-squares algorithm (see 
Sect. |3.2| ). The error on the abs olute line position is less than 
0.00094 cm“^ as derived in Sect. |3.3| The error on the FWHM 
is dominated by the uncertainty on the instrumental line profile 
as described in Sect. 13.41 The maximum error on the relative 
intensities due to uncertainties in the fiat-fielding procedure is 
estimated to be around 10% (see Sect. 1^. 

Column 14 contains information on problems, which were 
noticed during the data reduction process and could have an in¬ 
fluence on the determined line parameters. The information is 
given in the form of up to five letters that represent different cases 
as explained in Table[^ We call these letters fiags throughout this 
article. 

The number of spectral lines in the three spectra are summa¬ 
rized in Table|5 A total of 40408, 51776, and 58274 lines were 
detected and fitted for MW powers of 25 W, 50W, and lOOW, 
respectively. This confirms the trend we have already seen in 
Fig. the total number of detected lines grows with increasing 
MW power. Table also gives the number of identified lines for 
the different species. More lines from the literature are identified 
when the discharge is operated with a higher MW power because 
both the intensity of the lines increases (more fainter lines are in¬ 
cluded) and the lines become broader (makes a match with line 
positions from the literature more likely). 

The majority of spectral lines remains unidentified. Spectral 
modeling and calculated transition wavenumbers are needed to 
confirm that the unidentified lines in fact originate from nitrogen. 
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Table 4. Number of spectral lines in the three spectra. 



25 W 

50 W 

lOOW 

Total number of lines 

40408 

51776 

58274 

thereof lines without flags 

24585 

24863 

22723 

identified N I lines^ 

25 

27 

29 

thereof lines without flags 

17 

14 

17 

identified N 2 lines^ 

860 

953 

1003 

thereof lines without flags 

632 

636 

611 

identified lines^ 

886 

898 

903 

thereof lines without flags 

542 

476 

395 

unidentified lines 

38637 

49898 

56339 

thereof lines without flags 

23394 

23737 

21700 


Notes. Out of 29 lines from the NIST database. Out of 1260 lines 
of the (0 - 0) band of the B^IIg-A^Sy system of ^"^N 2 from 
( 1919\ . Out of 1017 lines of the (0 - 0), (0 - 1), (1 
bands of the A^Ifu-X^E^ system of ^"^N; 


from 


See main text and Table^for the meaning of flags 


Ferguson et al. 


Effantin et al. 
(if and (1-2 ) 


(19921. 


Because of the complex structure of the N 2 spectrum in the near- 
infrared, the theoretical modeling of this spectrum is beyond the 
scope of this paper. 

While the total number of detected lines increases with 
higher MW power, the number of lines that show no noticeable 
problems during the analysis decreases. This is mainly due to the 
many blends contributing to around 90% of the flagged lines (of¬ 
ten in combination with other flags). Total percentage of blended 
lines in the line list amounts to 35%, 47%, and 55% for 25 W, 
50W, and lOOW, respectively. 

4.2. Line intensities 

We present the distribution of line intensities using the spec¬ 
trum of the discharge operated at MW power of 50W. Figure 
displays the flat-field corrected line intensities Ic as a function 
of wavenumber. The corresponding raw spectrum is shown in 
Fig.|^ At lower intensities, the distribution exhibits an envelope 
with a shape related to the applied flat-field curve (compare with 
Fig. 1^. This means that our measurements are less sensitive to 
fainter lines toward larger wavenumbers, owing to the response 
curve of the InSb detector. It can therefore be expected that all 
emission lines with a peak intensity Iq > 0.4 a.u. are included in 
our spectra and line lists, but that the sample is incomplete for 
fainter lines. 

A histogram of the peak intensities 4 is shown in Fig. 
again exemplary for the discharge operated at MW power of 
50 W. The distribution peaks at around 0.1 a.u. and decreases 
quickly toward smaller intensified The shape of the distribution 
at intensities below 0.4 is biased by the wavenumber-dependent 
detection efficiency. The total dynamic range of the spectral lines 
covers about four orders of magnitude. Figure also highlights 
the fraction of lines with and without flags relative to the to¬ 
tal number of detected lines as a function of intensity. For all 
three measurements, the percentage of flagged lines stays above 
40% for the lower peak intensities between 0.01 and 0.1, and 
decreases to a value below 25% for intensity values above 10. 

Near-infrared detectors used for astrophysical spectrometers 
can usually distinguish line intensities over a dynamic range 
of two orders of magnitude. Therefore, only a fraction of the 
spectral lines of the molecular nitrogen spectrum are usable for 

^ The designation for arbitrary units, a.u., for the intensity is dropped 
in the remainder of the article for better readability. 
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Table 3. Explanations for the line list columns. 


Column 

SymboF 

Unit 

Explanation 

1 

2 

p 

o 

cm“^ 

designation of atom or molecule; “NoID” if not identified 
line center, calibrated as described in Sect. 3.3 

3 

Vo 

cm“^ 

line center 

4 

O’(vo) 

cm“^ 

error on vo, derived from weighted line fitting (“0.” if fit parameter is at boundary condition) 

5 

h 

a.u. 

line intensity, with fiat-field correction as described in Sect. 3.5 

6 

I 

a.u. 

line intensity 

7 

cr(I) 

a.u. 

error on /, derived from weighted line fitting (“0.” if fit parameter is at boundary condition) 

8 

O-Q 

cm“^ 

width of the Gaussian line profile 

9 

cricTc) 

cm“^ 

error on ctq, derived from weighted line fitting (“0.” if fit parameter is at boundary condition) 

10 

FWHMc 

cm“^ 

full width at half maximum, corrected for instrumental broadening as described in Sect. 3.4 

11 

FWHM 

cm“^ 

full width at half maximum 

12 

o-(FWHM) 

cm“^ 

error on FWHM, derived by error propagation from (t((Tg) 

13 


cm“^ 

distance of vq, c to reference line in literature (“99.” for not identified lines) 

14 

flags 


“x”: fit has failed or a fit parameter is at boundary condition. 

15 



“y”: fit residuals exceed at least 11 times (= intensity threshold-rnoise value) the noise 
level and are larger than 5% of the peak intensity. 

“b”: the neighboring line is closer than the sum of their FWHM and has more than 
half of the line’s peak intensity I (blend). 

“r”: FWHM is narrower than the resolution set for the FTS. 

“f”: (t(FWHM) is more than five time larger than the average error on FWHM (for 
unidentified lines) or FWHMc scatters more than five times from a linear trend for 
its species (for identified lines). 

transition as identified from Effantin et al. (1979) or Ferguson et al. (1992). 


Notes. as used in equations and text in this article. 


wavelength calibration of astrophysical spectrographs in prac¬ 
tice. However, as the utilization of a spectral line of molecular 
nitrogen for wavelength calibration should not depend on its in¬ 
tensity, an observer is flexible in selecting the appropriate expo¬ 
sure time for a certain calibration frame. 

4.3. Line density 

The spectrum of the nitrogen gas discharge exhibits emission 
lines over the complete spectral range under investigation. A his¬ 
togram of the line distribution is shown in Fig. for all three 
measurements. Around 350 to 1300 lines are located within one 
bin covering 100 cm“^. The number of lines in each bin increases 
with higher operating MW power. 

In order to be useful for wavelength calibration, a spectrum 
needs a sufficient number of lines, depending on the instrument’s 
resolution and wavelength coverage. Similar to Fig. 0Fig-0 
shows the distribution of lines as a function of wavenumber, but 
this time each bin covers a constant number of 1000 resolution 
elements with a resolution of R = /l/A/l = 100000. This corre¬ 
sponds to specifications of high-resolution grating spectrographs 
like, e.g., CRIRES ( |Kaufl et al.|2004| ). Considering the expected 
dynamic range of a near-infrared detector, in Fig. [T^ we only in¬ 
clude emission lines with peak intensities within two orders of 
magnitude: the histogram in the upper panel of Fig, contains 
lines with fiat-field corrected intensities of 0.5 - 50 (compare 
horizontal lines in Fig. |^. For example, at wavenumbers above 
9000 cm“\ an average of about 0.4 lines fall into one resolution 
element. This number decreases toward smaller wavenumbers 
until very few lines of this intensity are available at 4500 cm“^ 
To detect more lines at smaller wavenumbers, we need to go 
to lower intensities as shown in the lower panel, which can be 
achieved using longer exposure times and a longpass filter. The 
second histogram contains lines with intensities 0.05 - 5. Here, 


on average, there are approximately 0.5 lines per resolution ele¬ 
ment in the wavenumber range 5500-7750 cm“^. At wavenum¬ 
bers above 7750 cm“^ (dashed line) fainter lines are missing in 
our measurements. 

The continuous spread of more than 40000 detected lines 
over the wavenumber range 4500 - llOOOcm"^ ensures a suffi¬ 
cient number of lines for wavelength calibration. The exposure 
time might need to be adjusted when observing in different spec¬ 
tral regions to obtain the optimal number of lines, depending on 
a spectrograph’s specific design. 


4.4. The influence of different MW powers on line width and 
intensity 


Varying the microwave power supplied to the cavity changes 
the mean power absorbed by electrons from the electric field 
and therefore the rate of collisional excitations in the plasma 
( Jankowski & Reszke||2011[ chapt. 1.2). Here we present how 
the widths of the spectral lines become broader with higher MW 
power and how the intensities of the atomic nitrogen lines in¬ 
crease with MW power. 

Doppler broadening is expected to be the dominating effect 
on the line width after correction for instrumental line broading 
because of the low gas pressure and the Gaussian-shaped spec¬ 
tral lines. The theoretical description of Doppler broadening can 
be expressed as 


FWHMc 


8krin(2)^ 


-Vo, 


(3) 


where k is the Boltzmann constant, m is the mass of the species, 
and c is the speed of light (e.g., |Bernath||2005] chap. 1.3). Ac¬ 
cording to Eq. lighter species create broader spectral lines 
compared to heavier species (at a constant temperature T), and 
the line width increases linearly with wavenumber vq. 
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Table 5. Values of the temperature parameter T in Eq.j^as derived by 
linear regression. 



NI 

N2^ 


r(25W) [K] 

957 + 10 

681 ± 11 

809 ± 36 

r(50W) [K] 

1200 + 2 

825 ± 8 

967 ± 27 

r(ioow) [K] 

1339 + 2 

970 + 7 

1113 + 20 


Notes. 1 band of the B^rtg-A^S^ system. 4 bands of the A^IIu- 
system. 


The FWHMc of all lines without a flag is plotted against 
wavenumber in Fig. [m The identifled lines of atomic and 
molecular nitrogen are highlighted. A linear regression using 
Eq. is applied to the data points of each identifled species, 
as represented by the lines in the respective colors in Fig. [TT] 
The linear regression weights each point according to its un¬ 
certainty in FWHMc, which can lead to a slight displacement 
of the linear regression relative to the bulk of the data. Atomic 
nitrogen clearly follows a trend of larger FWHMc as com¬ 
pared to molecular nitrogen. We also see the expected linear 
trend with wavenumber for all species. The distribution of the 
FWHMc from the unidentifled lines is consistent with a com¬ 
bined FWHMc distribution of lines from N 2 and NJ with a ratio 
of 0.82 : 0.18, respectively. The overall scatter is, however, too 
large to make conclusive statements about the origins of individ¬ 
ual lines. 

Equation includes a temperature parameter T. The defi¬ 
nition of a single gas temperature is here not possible because 
a microwave induced plasma is not in a local thermodynamic 
equilibrium. Each species in the plasma rather has its own tem¬ 
perature ( [Jankowski & Reszke||2011| chapter 1.4). Because a 
change in temperature for a certain species is related to a change 
in FWHM, we compile the values of T as derived from the lin¬ 
ear regression for completeness in Tablej^ The temperature rises 
for all species if the MW power is increased. The error on T, as 
obtained directly from the fit, is very small (< 1 K) because the 
linear regression consists of only one free parameter. To esti¬ 
mate more realistic uncertainties, we varied the data points by 
+(t(FWHMc) and report the resulting change in T in Table 

Furthermore, we analyzed the change in peak intensities for 
the 25 transitions of N I that are identifled in all three spectra. 
Figure displays the mean change in peak intensity when the 
MW power is increased from 25 W to SOW and then to lOOW. 
The values are calculated in the following way: the intensity 4, p 
of each NI line in the spectrum of the discharge operated at MW 
power P is divided by the intensity /c, 25 w of the same line in 
the spectrum of the discharge operated at 25 W; then, the mean 
value is calculated. The latter appears to follow a linear trend 
as indicated by the dashed line in Fig. and described by the 
following linear regression: 


IAlL.\ = 0.0756 — • P[W] - 0.8679. (4) 

\h,25W I W 

The individual intensity values are well separated into two 
disjointed distributions for the measurements with 50 W and 
lOOW, with only one outlier at 7349 cm“^. If we remove this 
outlier, the error bars shrink by about a factor of two and the 
linear regression becomes {p^) = 0.0781 ^ • P[W] - 0.9318. 

Article number, page 8 ofpfl] 


5. Conclusion 

We have presented the characterization of three emission spec¬ 
tra of a nitrogen gas discharge in the range 4500 - llOOOcm"^ 
created by a microwave-induced plasma. A speciflc character¬ 
istic of the experimental setup is the sealed gas cell, which al¬ 
lows one to repeat the measurement with the same gas proper¬ 
ties. The spectral measurements were carried out for three dif¬ 
ferent MW powers (25 W, 50W, and lOOW). We detected and 
fitted 40408, 51776, and 58274 emission lines in these spec¬ 
tra, respectively. Spectral line lists with line parameters are pro¬ 
vided for all detected lines. The spectra exhibit densely-spaced 
lines throughout the whole wavenumber range under investiga¬ 
tion with about 350 - 1300 lines per 100 cm“^ Depending on 
the MW power, between 35% and 55% of all lines are blended. 
We identifled molecular and atomic nitrogen lines using existing 
line lists from the literature. The majority of the lines remains 
unidentifled. Theoretical transition wavenumbers are needed to 
confirm whether they originate from nitrogen and to properly as¬ 
sign all transitions. 

The analysis of line intensities and line density focused 
on the potential application as a wavelength calibration source 
for astrophysical spectrographs. The number and distribution 
of emission lines qualifies the spectrum of the recorded source 
to be used for wavelength calibration over the whole spectral 
range under investigation. For a resolution ofR= 100000, ap¬ 
proximately one line every two resolution elements can be real¬ 
ized over a spectral range of a few thousand wavenumbers with 
peak intensities within two orders of magnitude. Depending on 
the wavelength range of a speciflc spectrometer, exposure times 
need to be adjusted to select, in combination with a bandpass 
Alter, a suitable subset of emission lines. Tests for stability and 
aging effects need to be carried out to further investigate the fea¬ 
sibility of the molecular emission spectrum as a reliable wave¬ 
length calibrator. 
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Fig. 6. The recorded spectra between 9892 and 9895 cm ^ From top 
to bottom: discharge operated at MW power of 25, 50, and 100 W. The 
blue line is the fit composed of a superposition of Gauss profiles to the 
peaks exceeding the minimum peak intensity threshold (dashed hori¬ 
zontal line). A green square indicates an identified N I line, while a 
gray circle marks an unidentified line. See main text for more informa¬ 
tion (Secf |3.8| ). Note that the ordinate axis is not the same scale for each 
plot. 
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10 ^ 



wavenumber [cm“^] 

Fig. 7. Line intensities after flat-field correction as a function of 
wavenumber. Note the logarithmic scale. The horizontal lines indicate 
intensity ranges used in the upper panel (dashed lines) and lower panel 
(dash-dotted lines) of Fig. Data from discharge operated at MW 
power of 50 W. 



flat-field corrected line intensity [a.u.] 

Fig. 8. Histogram of the peak intensities Ic of the detected lines. Bars 
are stacked on top of each other. Data from discharge operated at MW 
power of 50 W. 



Fig. 9. Histograms showing the distribution of emission lines as func¬ 
tion of wavenumber. Each bin covers a range of 100 cm" ^ The bars are 
not stacked and show the absolute number of lines in one bin for each 
of the three spectra operated at different MW powers. 



wavenumber [cm“^] 


Fig. 10. Histograms of detected emission lines as a function of 
wavenumber. The samples are limited to emission lines with a peak 
intensity 0.5 - 50 (upper panel) and 0.05 - 5 (lower panel). Each bin 
covers 1000 resolution elements based on a resolution of R = 100000. 
In the lower panel, the dashed vertical line is located at 7750 cm"^ and 
the sample is expected to be incomplete at higher wavenumbers (com¬ 
pare with Fig.[^. Data from discharge operated at MW power of 50 W. 
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n: 
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■ N I 
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wavenumber [cm“^] 
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V (A-X) 


MW power = 50 W 


•.Ar 


0.12 
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0.08 

E 
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0.04 

0.02 
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Fig. 11. Each point represents the FWHMc of one spectral line as a 
function of wavenumber. The solid lines show linear regressions with 
ordinate crossing forced at the origin. The black dashed horizontal line 
indicates the FTS resolution. From top to bottom: discharge operated at 
microwave power of 25, 50, and 100 W. 


wavenumber [cm“^] 


• No ident. 
■ N I 




Fig. 12. The mean intensity of the detected N I lines relative to the 
respective line intensity in the spectrum of the discharge operated at 
25 W. The green dashed line indicates a linear regression. 
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